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Abstract 

The pre-equilibrium contribution to open charm production in nuclear col- 
lisions at yfs = 200 AGeV is calculated using three different models for the 
correlations between momentum and space-time coordinates. Ideal (Bjorken) 
correlation between the rapidity y and space-time rapidity r] of mini-jet gluons 
suppresses greatly the pre-equilibrium yield and even allowing for the min- 
imal uncertainty correlations leads, in contrast to previous estimates, only 
to a small pre-equilibrium charm yield as compared to initial yield due to 
gluon fusion. The "intrinsic" charm process is negligible in the mid-rapidity 
domain. 



I. INTRODUCTION 



Open charm production, direct photon, and dilepton production are among the most 
direct probes of the early time evolution of the quark-gluon plasma produced in ultra- 
relativistic nuclear reactions. At collider energies ^/s > 200 AGeV the initial mini-jet 
plasma is mostly gluonic |5],|6|] with a quark content far below its chemical equilibrium value. 
Furthermore, the initial transverse momentum distribution of those gluons is very broad |2j 
resembling a hot thermal gas of gluons with an effective temperature T ~ 500 MeV ||. 
Because charm is produced mainly through gluon fusion, open charm production provides a 
probe of that initial gluonic state. In contrast, hidden charm |7| is mostly sensitive to final 
state interactions in the later stages of evolution. Photons and dileptons are complementary 
probes of the evolution of the suppressed quark component of the plasma. 

The present study is motivated by two recent studies p],|4| of open charm which predicted 
widely different rates in nuclear collisions. In ref. [[| the pre-equilibrium contribution was 
found to be almost equal to initial gluon fusion rate. A similar factor of 2 enhancement of 
charm from thermal production in hot-glue scenario was also suggested ||. In ref. [|J, a 
more provocative claim was made that open charm may even be enhanced by over an order 
of magnitude above the initial pQCD rate. The main result of our present study is that 
correlations between the rapidity y and the space-time rapidity r\ lead to a large suppression 
(about a factor of 40) relative to the uncorrelated case. Thus, the pre-equilibrium open 
charm production is found to be unfortunately a very small fraction of the initial fusion 
rate. The large enhancement of charm production in ref. @ is found to be due to an 
overestimation of the contribution from the flavor excitation processes and the use of a low 
energy A a scaling from pp reactions measured at Ei a b = 300 — 400 GeV. 

The paper is organized as follows: In section |J the dependence of the direct pQCD 
rates for charm production on structure functions, Q 2 scale, and K factor is reviewed and 
compared to existing data. The beam energy dependence and the A dependence of the 
initial charm production are compared to results in ref. [[| . In section [TIT| , the pre-equilibrium 
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charm production is calculated. The mini-jet rapidity and transverse momentum distribution 
are fit to results of the Monte Carlo HIJING model || including initial and final state 
radiation. Three different models for the space-time and momentum correlations are studied 
and the influence on the charm yield is calcualted. Of the three models, we concentrate on a 
minimally correlated model resulting from uncertainty principle, which is similar to the type 
of correlation assumed in ref. ||. We also study the sensitivity of the results to different 
models of the formation physics . Section fV| contains the summary. 



II. INITIAL CHARM PRODUCTION 



Heavy quark production in pp reactions was calculated long ago in ref. ]TT[ including both 



fusion and heavy flavor excitation processes in the leading order pQCD. It was proposed that 
the flavor excitation processes were dominant at high energies because a small Q 2 exchange 
can easily liberate any charm component in the nucleon while gluon fusion was suppressed 
because Q 2 > AM 2 . In the Parton Cascade Model |Q, both mechanisms are incorporated to 
calculate s, c, b quark production in nuclear collisions. There the results suggested that the 
flavor excitation of the charm quark of nuclear structure functions would be the dominant 



source of charm production in nuclear collisions as well. However, it is pointed out [12 
that the original flavor excitation rates in ref. |ll]] were too high in the Xf ~ region 
due to neglected interference with other pQCD amplitudes to the same order. When all 
diagrams were added together, a large destructive interference was found to suppress the 
flavor excitation rates by powers of A/M q , where A ~ 300MeV is a typical QCD scale and 
M q is the heavy quark mass. The suppression factor appears to the process g + c(c) where the 
charm is evolved from the structure functions using perturbative QCD, as also shown in ref. 
[1"3| . We note that there is a possible non-perturbative charm component (intrinsic charm) 
in the nucleon. There are experimental constraints on the amount of that non-perturbative 
charm component |I^|T5|| . The total contribution of the intrinsic charm was shown in refs. 



TTj| , [T7[ ] to be small (about 10%) in the midrapity region where most of the charm is made. 
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Although the contribution of the intrinsic charm component appears important at large Xf, 
its contribution to the total cross section is small and well within the uncertainties from 
other sources. 

In this paper we only include fusion processes for the parton level cross sections as in 
ref . |2J . For the production in p — p collisions, we use the light quark and gluon structure 
functions from Gliick et al. fl8| and Duke-Owens |19| for comparison. The pQCD differential 
cross sections for a + b — > cc + X, are taken from ref. (T^j. For example, 
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where \ = y 1 — 4M 2 / s and we consider the following two choices for the scale Q 2 in the 
coupling constant a s (Q 2 ) = I2ix / [(33 - 2n f ) log(<3 2 /A. 2 )] from ref. [[TTJ] : 

1. for gg — > cc, Q 2 = s/2; for qq — >■ cc, Q 2 = s. (Q 2 choice-1) 

2. for both gg — >• cc & gg — > cc, Q 2 = s. (Q 2 choice-2) 

We take n/ = 4 for charm quark production and n/ = 5 for bottom quark production. 
The QCD scale A depends on the choice of parton distribution functions and is given in the 
tab 



e below 


Parton distribution functions 


A(GeV) 


GRV-LO set 


0.25 


GRV-HO set 


0.20 


Duke-Owens set l(DOl) 


0.20 


Duke-Owens set 2(D02) 


0.40 



To incorporate approximately the next-to-leading-order corrections to the above rates we 
multiply the leading order results by a K-factor. In general, K-factor depends on the choice 
of parton distribution functions, the center of mass energy of the collision, and the type of 
the projectile and target particles. Calculations to order 0(a^) for the subprocesses were 
carried out |20|j2~l[| , and afterwards the calculations to order 0(a^) foip + p collisions were 
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made [22,23]. For D01, M c = 1.5 GeV, Q 2 = 4:M 2 ,P lab = 100 - 1000 GeV, the K-factor 



for p — p collisions was found to range from 2.85 to 4.1. We also note a recent result 



24[ where the dependence of the K-factor on the final momentum of the initially produced 
charm was studied for high energy AA collisions. As a function of the rapidity of the charm, 
the K-factor is almost a constant ~ 2. As a function of p±, the K-factor increases from 1.3 
at p ± = 0.7 GeV to 3.4 at p ± = 6 GeV. 

In Fig. 1 we compare the so calculated charm cross section to the limited data on inclusive 



cc production in p — p collisions. The NA34 data for a c h ar m is taken directly from ref. |25 
The values for the other data lines are computed from D-meson cross sections according the 
argument in ref. [^6| by using the published experiment results [p?|-|30[|. Earlier experiment 



results |3T| also show big uncertainties among the different experiments. 

In Fig. 1 We see that both the solid curve and the dashed curve fit the low energy 
data reasonably well, so we use these two parametrizations for the following high energy 
calculations in this section. As a consistency check, we also plot the long-dashed curve using 
the same parameters as in Fig. 1 of ref. |£2J (i.e. DOl, M c = 1.5GeV, Q 2 = 4M 2 ) using 
constant K = 3 for simplicity. Comparing the solid and dot-dashed curves shows the strong 
dependence on the assumed charm quark mass for the GRV-HO set. Comparing the solid 
and dashed curve we see that different choices for the Q 2 scale can be compensated for by 
shifts in the K factor. These results together with the large uncertainty of data emphasize 
the need to measure pp and pA to fix uncertainties in the initial charm production rate in 
order that charm production in AA can be properly calculated. 

Next we compare our results for the rapidity density of produced cc pairs at Y = with 
results of ref. [[|. In Fig. 2 the energy dependence in the range between RHIC and LHC 
( ^/s = 200 — 6300 AGeV) for Au + Au collisions are shown. The scaling from pp results to 
AA is 



where crfneiastic * s taken from ref. |32| . Glauber geometry for central high A + A collisions 



gives a = 1.. In Fig. 2 the solid curve is our result using the same parameters as for the solid 
curve in Fig. 1. The parametrization for the dashed curve in Fig. 1 gives a curve higher than 
the solid curve by 15% to 30%. The four long-dashed curves, curvel to curve4, are all from 
PCM calculations The top curve4 is the parton cascade model result for the so-called 
QGP formation case, including both the fusion and the flavor excitation processes. That 
curve is higher than our solid curve by about an order of magnitude because it includes the 
contribution from flavor excitation processes. Curve3, the curve with filled squares, shows 
the contribution to curve4 from fusion processes only (processes (1) and (2) in the notation 
of ref. Q), and curve3 is very close to our results. The bottom curvel is the estimate without 
QGP formation by extrapolating the parton model pp result to AA using A im scaling. It is 
lower than our solid curve by a factor of 6 to 2.5. The main source of this difference is from 



the A-dependence of p — A cross sections. Ref. [§]] used an A a scaling with a = 0.76 |33 
instead of the value a = 1 we use from Glauber geometry. We note that the value a = 0.76 is 
taken from low energy experiments, where energy conservation suppresses the contribution 
from multiple collisions. At high energies, QCD factorization implies that a = 1 for p — A 
scaling is the appropriate scaling modulo small nuclear shadowing effects. To demonstrate 
this effect from different A a scaling, we multiply Curvel by a factor of A 10 /A ' 76 = 3.55 
and get Curve2, which is close to our results. In summary, the factor ~ 50 enhancement of 
charm production suggested in ref. [|J comparing curvel with curve4 for charm production 
at RHIC is a consequence of the inclusion of incoherent flavor excitation processes and the 
extrapolation from pp to AA via low energy scaling. Given the coherent suppression of the 



flavor excitation processes |12| and the high energy scaling under consideration, it is only 
sensible to compare curve2 with curve3. In that case Fig. 2 leads to the expectation that 
the pre-equilibrium charm production should be comparable to the initial fusion rate. This 
removes the bulk of the discrepancy between ref. and ref. [f|]. 

As a further check on the parameters we compare charmed hadron Xf results in Fig. 3 



with 400 GeV p — p data [27J using the idealized 5-functionfragmentation function. The 



realistic fragmentation function used in ref. |L6[ lowers the curves slightly and reveals the 



true high-x/ intrinsic charm component. In Fig. 4 we compare bb production. Here we take 
M c = 4.75 GeV as in ref. J23|, with K = 3,n f = 5. The data point at y/S = 630GeV is 



from ref. [33J]: a(pp — > b + X) = 19.3±7(exp.) ±9(th.)fj,b, and only the experimental error is 
indicated in Fig. 4. At Vo = 1.8Tet>, our value is 41.8/ife x K = 125/ife. This is significantly 



larger than found in ref. [23] 



III. P RE-EQUILIBRIUM CHARM PRODUCTION 

We consider next the pre-equilibrium contribution to the charm yield in A + A. This is 
the charm produced through final state interactions between partons in the dense mini-jet 
plasma. Here we only calculate the dominant contribution from mini-jet gluon fusion. 

A. Spectrum of Mini-Jets 



The spectrum of mini-jet gluons in leading-order follows from ref. |35| 
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The term mini-jets refers to unresolved jets at a scale p± > p± cut 
cross section to produce mini-jets is given by 



(4) 
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2 GeV. The inclusive 
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dydp\ 

where f\ is the incident parton distribution evaluated at x± = p±(e y + e m )/y/s at a 
scale Q 2 = p± 2 . The light-cone coordinates of the initial and final partons are p\ = 



2xip ,0,0 



P2 



0,2a; 2 po,0 , p 3 = [m_|_e y3 , m±e V3 ,—p~±], and the observed parton has 



P = [m±e y ,m±e y ,p±\. The subprocess Mandelstam variables are s = SX1X2 etc.. For the 



calculation of mini-jet gluon fusion process in the following section |111 BL we choose Q 2 = s. 



As in ref. 0, we use DOl as the proton structure functions and K = 2, M c = 1.5GeV 



for the mini-jet production. Shadowing on Au is taken from ref. p6[. The resulting trans- 
verse momentum distribution of mid-rapidity mini-jet gluons at ^Js = 200 AGeV is shown 
by the open circles in Fig. 5. We call this distribution the hard distribution since it has 
Vi-cut = 2GeV. It is compared to the solid line, which is the output of the Monte Carlo 
calculation via the HIJING model || that includes initial and final state radiation. 
For convenience we have parameterized the Monte Carlo results as the following: 

dN 



g(p ± )p(y)A A/3 = 0.06e- L25p ^ cos 



A 4/3 , with \y\ < 3.7 (7) 



dydp± 

In the following, we call this parameterized distribution the soft+hard distribution. The 
soft+hard, hard, and Monte Carlo distributions are very close to each other in the semi- 
hard p± > 2 GeV region at y — 0, as seen in Fig. 5. However the parametrized distribution 
falls underneath the Monte Carlo result in the region p± < 1 GeV. We emphasize that the 
soft component is strongly model dependent as it requires the furthest extrapolation from 
the pQCD hard domain. The Hijing yield in that region is due to initial and final state 
radiation. Other contributions in this soft domain from coherent string are possible ||. 
While most of the following results are obtained with the simple parametrization above, we 
will check the sensitivity to variations of the soft component as well. We also note that 
at bigger rapidity the p± spectrum falls more rapidly. The above parametrization does not 
include that property. However, that property only lowers the high p± tail, and hardly 
changes the low p±_ part and the total number of the pre-equilibrium charm. 



B. rj - y correlations 

1. Bjorken correlation 

In ideal Bjorken dynamics, the space-time rapidity rj = 1/2 log[(t+z) / (t— z)\ and the true 
momentum rapidity y = l/2\og[(E + p z )/(E — p z )] are assumed to be perfectly correlated. 
This is referred to as the inside-outside picture and the phase-space distribution function in 
this case has the form 
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f <*-« *>« = Mi£k s( * - y)e(T - n)e{Tf - t) (8) 

Ti = 0.1 fm/c is the mini-jet formation time. Tf pa 1.7 fm/c is the proper time when 
the energy density of the pre-equilibrium mini-jets falls by an order of magnitude to ~ 
2GeV / fm 3 due to rapid longitudinal expansion, and that is when we terminate the pre- 
equilibrium stage. 

The phase space distribution is normalized such that 

r F(x,p,t) B jd 3 x _ d 3 N _ 1 dN 

J {2n) 3 ~ ~i?p~ ~ E dydpl ^ 
In this section we study the pre-equilibrium charm production at y — [§]: 

( B wL -Mm^^^ m |M|2 m 

Denoting dN/dydp± = g(y,p±), and p± 1 = (cos0i, sin0i, 0)pxi, the ideal r]—y correlation 
leads to 



d 3 pj y=0 4 32(2nfT7cR\J ' ^ ^ ^ ^ ^ 



32(27r) 2 7r_R^ 7 p_i_ 2 [1 — cos(0i — 02)] — (E cosh?7 — pcos^i 

In deriving the above, we have used kinematic relations 

E' = + p± 2 ) coshr? — E 
6&E) S(p±i-P±i,o) 



E P±2 [1 ~ COs(0i — 02)] — (-EcOsllT? — pCOS0i) 

p ±2 (E coshr] -pcos0 2 ) ^ 

,0 p± 2 [1 — cos(0i — 2 )] — (E coshr/ — p cos0i) 

Numerical integration of the above integral in equation (|Tlj) leads to the results shown in 
Fig. 6. The solid line is the p_|_-distribution for the initial charm production, from section [II]. 
We see that the pre-equilibrium contribution in this strongly correlated case is totally neg- 
ligible. This result is similar to the thermal charm production contribution calculated in 
ref. except that in our case the curve extends to higher p± because of the broader initial 
mini-jet distribution in p±. 
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2. Uncorrelated r] — y 



In ref. [0], another extreme case, opposite to the ideal Bjorken picture, was considered. 
In that case the gluon distribution is assumed to be completely uncorrelated as in an ideal 
thermal fireball. This assumption leads to 

r,/^ - X ( 27r ) 3 1 <IN 

Fmt)n= LJ--— (13) 
If one assumes a fixed volume V = t(kR a , then / dt ~ -ry — Tj, and 

^ - 1 T/ as in ref. [2]. (14) 



V 2 ttR 2 a Ti 
Then from equation ([K]), we have 
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Hp 



dy 1 dy 2 dp ±2 d(j) 1 d(j)2 g(yi,P±i,o)g(y2,P± 2 )\M\ 2 

cosliyi cosh?/2 P_i_2 [cosh(yi — 1/2) — cos(0i — ^2)] — (-E cosIm/i — pcos^i) 

^(p±i-p±i,o) 

p_L 2 [cosh(|/i - 1/2) - cos(0i - 2 )] - (E coshyi - p COS^i) 

p± 2 (E coshy 2 -pcos0 2 ) , . 

11,0 P±2 [cosh(2/i - y 2 ) - cos(0i - 2 )] - (E cosliyi - p cos</>i) 

For uncorrelated case, the pre-equilibrium charm production is much larger than the Bjorken- 

correlation case, and is comparable with the initial charm yield, as shown in Fig. 7. This is 

similar to the result in ref. where the pre-equilibrium charm production has almost the 

same magnitude and p^-shape as the initial charm. 



3. Minimally- correlated rj — y 

We consider here the simplest source of i] — y correlations resulting from the minimal 
geometrical spread in initial production points required by the uncertainty principle. This 
type of correlations are included in the parton cascade model and discussed in ref. JJ. The 
phase space distribution function including such minimal correlations has the form 
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F(x,p,t) M in = A/" / ^\ ' ~~7777^~ po{x Q ,t )6{x - x - vAt)d 3 x dt . (17) 
y dydpi i + (Jgi)2 

The integration is over the space-time coordinates (a?o, to) of the production points of the 
gluons. These points are distributed according to a normalized density p (x ,t ). The delta 
function arises to take into account the free streaming of the partons from the production 
point, with velocity v = p/E, where E = p± coshy and p z = p± sinhy. The theta function 
defines what we mean by pre-equilibrium. The proper time when the pre-equilibrium fusion 
is terminated is T max , which is determined below in Fig. 8. The theta function insures that 
only those gluons with proper time less than r max contribute. 

The formation physics is included via the Lorentzian formation factor fH| 



[i+(t / (p)/Af)r i , (is) 

where At = t — t is the elapsed time, and the formation time is given by 

, 0.2GeV, r , 

tf(p) — cosm/ \J m ) ■ (19) 

P± 

We note that the above formation factor more accurately describes the interference phenom- 
ena suppressing production at early time than the conventionally assumed factor 

9[At-t f (p)} . (20) 

In the following we consider both formation functions for comparison to check for the sen- 
sitivity to this formation physics. 

We assume that / po(x~o,to)d 3 xodto = 1. In this case the normalization factor is M = 
(2n) 3 /E, so that 

lim f F{x,p,t) M ind 3 x/{27r) 3 = d 3 N/d 3 p . (21) 

As discussed in ref. 0, the production points are spread along the beam axis according to the 
uncertainty principle by an amount 5z = d ~ fr/p± since the dominant parton interaction 
leading to a y = parton with final p± has an initial longitudinal momentum xPq ~ p±. 
We take particular model 
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d=^(/m) (22) 

Clearly this is only a rough guess, but it allows us at least to investigate the sensitivity of 
the results to a particular r\ — y correlation that results from this spatial spreading of the 
production points. We emphasize that it is precisely the uncertainty of the initial space- 
time formation physics that leads us to study the possibility of open charm production as 
an experimental probe of that physics. 
Given the above assumption we take 

Po(xb,t Q ) = — B2-£(*o) — t=t~ (23) 



where d is the mean spread for gluons depending on p± from above. This distribution only 

spreads out the production points along the beam axis. A more realistic treatment would 

also smear out in the time coordinate. 

Neglecting transverse expansion, we obtain finally 

F(¥ * ,\ _ (2?r) 3 P± (z _ tanhyt?( p± f/2 1 dN 0{r max -^—) 
[,P,)Mm ~V^7rR\0.2 e pdydpll + C-^y)^ [M) 

Let ai = tanhyi, a 2 = tanhy 2 , &i = (%^) 2 /2, b 2 = (^f) 2 /2, then after integration over z, 
we have the final expression as the following, while its numerical results are shown in Fig. 9: 



d 3 P J y= o IQ(2it) a R\ J p_L 2 cosh?/! cosh?/2(£'cosh?/2 — pcos0 2 
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0.2 cosht/2 *\ 2 



1+ v 



(25) 



In the above tf = T max min (coshyi, coshy 2 ), and p± 1 is the same as in equation fll6|) . Note 
that by using the unit GeV^ for momentum and unit fm for time, the expression 
in equations ([11]), (0), and ( f25|) has the dimension GeV~ 4 fm~ 2 , and we need a factor 
(he) 2 ~ (0.2GeK/m) 2 to convert it to the dimension GeK -2 , which we have used in Fig. 6, 
Fig. 7,Fig. 9 and Fig. 10. 

We also plot the energy density curve at z = as a function of time in Fig. 8. We see 
that it increases first, and reaches maximum at the time about 0.1 fm/c, then the energy 
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density decreases linearly to ~ 2GeV/ fm 3 at ~ 0.9/m/c ( 1.7 fm/c ) for hard (soft+hard) 
distribution. We choose the above time as the cutoff r max . 

The previous uncorrelated case neglects the finite formation times of the mini-jets. In 
order to see the formation-time effect, we also use the ^-function form in equation (EDI) 



instead of the Lorentzian form in equation (TT8|) for the formation-time effect. The result 



from this ^-function is about 10% higher at p± = OGeV, and 10% lower at p± = 9GeV, as 
shown in Fig. 10. The lack of sensitivity to the formation-time physics is due to the relative 
large p± for the gluon mini-jets in the charm production process. There would be more 
sensitivity had the production been dominated by low p± components. 

We also see that for the soft+hard distribution the soft gluons significantly increase 
the pre-equilibrium charm production in both low-p^ and high-p^ region, with the largest 
increase in low-p^ region. It is interesting to identify where the enhancement comes from. 
In Fig. 9, the curve with diamonds shows the contribution from the fusion of soft gluons 
both with p^ < 2 GeV, and the curve with unfilled squares shows the contribution from the 
fusion of hard gluons both with p±_ > 2 GeV. These two curves are both very low compared 
with the curve calculated from the soft+hard distribution. So the enhancement going from 
hard distribution to soft+hard distribution mainly comes from the fusion of hard and soft 
mini-jet gluons. 

We have noted before that our fit for the mini-jet gluon spectrum falls below the Monte 
Carlo result from HIJING calculation. We can fit the soft gluons from HIJING better 
by using 0.265e~ 2 6pj - for p ± e (0, 1.1) GeV, and use the old fit 0.06e~ L25pj - for higher-p^ 
gluons. This new fit gives us more very soft gluons. We have done the calculation for 
minimally-correlated case using the new fit, and the result is different only by less than 
10%, which means the super-soft gluons are not very important for the pre-equilibrium 
charm production. 

There is also a possible cross-term contribution from the interactions of the incoming 
nuclei and the pre-equilibrium gluon mini-jets. However, our preliminary result shows that 
it is not larger than the above pre-equilibrium charm yield and is therefore also negligible 
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compared with the initial charm production. 



C. Why is the pre-equilibrium charm yield so small? 

To understand the reason why the pre-equilibrium charm yield is so small compared to 
the initial yield as found through tedious numerical calculations in the previous section, 
we consider here the calculation of the total number of pre-equilibrium charm pairs. The 
expression for that number is given by 

N = 4^ I * X I ^^nx,Pi,t)F(x,p 2 ,t)sa(s) (26) 

where a(s) is the integrated cross section for the process gg — > cc, see equation(H). Our main 
strategy is to estimate the mean difference between the two gluon rapidities, then from the 
kinematical constraint on charm production (s > 4:M^) estimate the effective lower cutoff 
for p^ of the mini-jet gluons. Thus we separate the p± integrals from the rapidity integrals, 
and have a rough estimate for the total number of charm pairs. 



1. Bjorken correlation case and Uncorrelated case 
For the fireball case, 

, (2tt) 3 1 dN 
F(x,p, t) Fb - 



p V dydp~± 

s = 2p ±1 p ±2 [cosh(yi - y 2 ) - cos(0i - </» 2 )] (27) 



For the Bjorken case, 

- \ (2vr) 3 dN SOn-y)^, 
F(x,p, t) Bj = y J W f 9r - n 6 (r f - t) 

p± dydp± Tirri A 

s = 2p ±1 pj_ 2 [1 - cos(0i - 4> 2 )} (28) 

For all the cases, we use the fit to the gluon distribution given by equation(^), where g(p±) = 
ae -b P± = .06e- L25pi . Therefore, 
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Npb = " IttR 2 / dyi coshyi / dV2 m^m I dp±l9 ^ P±1 ^ I dp±29 ^ p± ^ J <tyi J d<f> 2 sa(s) 

(hcfln^A 8 / 3 f f f f f 

N B j = 1^ J drj [p(r])} J dp ±1 g(p ±1 ) J dp ±2 g(p ±2 ) J dfa J d<f) 2 sa(s) (29) 



The dominant contribution is coming from the vicinity of the production threshold where 
s = 4M 2 = 9GeV 2 fTT|, so we make the following rough estimates: 



sa(s) ~ a (s) ~ 0.06 
y dp ± g{p ± ) ~ y^ dp ± g{p ± ) ~ (30) 

where p c is the effective cutoff value for p± 1 and p± 2 from the requirement s > 4M 2 . 
For the fireball case, 

Jd yi dy 2 ^±^2^ cosh(y 1 - y 2 ) 
(<»*<*-*)) = f dyidy2 ^L^L 4 '° < 31 > 

j J^cosnyi cosnj/2 

Since the mini-jet p± spectrum is dropping almost exponentially, the production heavily 
favors the smaller cutoff p c , so the mean value of cos(0i — <p 2 ) is most likely to be negative. 
We take (cos(0i — 2 )) ~ —0.5. Then s ~ 9p 2 , so the effective cutoff for the fireball case is 
p c ~ 1.0 GeV. 

On the other hand, for the Bjorken case, 

Vi = V2 = V ^s~3p 2 c =*p c ~1.73GeV (32) 



Using the same values as in section [III B| : T{ = 0.1 fm, Tf = 1.0 fm for fireball case, 
1.7 fm for Bjorken case, and 

p(yi) 



h^r 2 - 8 - /*'W")i 2 ~« < 33 > 



We then have the estimate for the total number of the pre-equilibrium charm: 



N Fb ~ ~^r- 2-8 2 (^) 2 (2vr) 2 [sa(s)) ~ 3.5 

iV fli ~ 4.9(^;) 2 (27r) 2 [$*(*)] ~ 0.098 (34) 
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Therefore we estimate Nph/^Bj ~ 35, in rough agreement with the detailed numerics. We 
see that the main source of the large increase going from the Bjorken case to the fireball case 
comes from the different p± cutoff. In the uncorrelated fireball case, one allows particles 
with different rapidities to interact with each other (see equation (|3~ID), thus more low p± 
gluons can take part in the interaction. Since the mini-jet p± spectrum is dropping almost 
exponentially, the fireball case produces a lot more pre-equilibrium charm than the Bjorken 
case (a factor of 6 increase from the smaller p± cutoff). Although the questionable linear 
proper time dependence in the fireball case also gives an considerable increase (about a 
factor of 3.5), it is not as important as the correlation effect. 



2. Minimal correlation case 

For the minimal correlation case, the estimate is unfortunately not as straightforward. 
The phase space distribution function is 

t?(--+\ (2tt) 3 dN t nccoshy 

F(x,p,t) M in = r— „ 9 Z T , , &{l~max 7~ )&{ T max ) (35) 

V27r7ri?j, nccoshy ay dp± cosh?/ p± 



and s is the same as in equation(|2T|). In the above distribution function we choose to use 
the ^-function for the formation-time effect. We have seen from Fig. 10 that Lorentzian 
formation-time formula and ^-function formula give almost the same result. 
Using equation (pp]) and after the integration over z, we have 

(he) 2 A 8 / 3 



[Tic) A ' f f f 

N M in = 47rR2 r- Jdp ±1 g{p ±1 ) J dp ±2 g(p ±2 ) J dfa Jc 



xjdyi — : dy 2 — - — sa{s) / dt = (36) 



coshy 1 J cosh?/2 Jt min yjl/h + l/b 2 
where 

t m in = hemax ( C ° Sh Vl , V2 \ , t max = T f min(coshy 1 , coshy 2 ) (37) 

V P±l P±2 J 

and ai,a 2 ,bi, b 2 are defined the same as in equation (|2"1|). 
We estimate that for the dominant part of the integral 
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11/ 2hc\ he coshy 



] ' i i i ,*iron , t max ~ 7"/ COsh y (38) 

01 &2 V / Pc 

where y = (\yi\ + 1 2/2 1 ) /2 . Now let u — t p c / (he coshy), then the last 3- dimensional integral 
in equation(|36D without the factor sa(s) is 

J ~ fdv,-^- [dy 2 -^ 1 -coshy du e - ^ ' V (39) 



1 cosh yi J 2 cosh y 2 2 ji 

The w-integral gives a dependence on Tf which is similar to the logarithmic dependence in 
the Bjorken case, and the exponential form in the integrand forces the spread y\ — y 2 to be 
small. Numerically, by taking 77 ~ 1.7 fm/c , p c ~ 2.0 GeV (as the first-step value) in the 
■u-integral, the above 3-dimensional integral is J ~ 19.1 , and when the integrand is weighed 
by cosh^x — y 2 ), the integral is ~ 23.3. So 

(cosh(yi - y 2 )) ~ 23.3/19.1 ~ 1.22 s ~ 3.Up 2 c => p c ~ 1.62GeV (40) 

Note that the above determined value of p c is insensitive to the first-step p c value we tried 
in the u-integral. 

Therefore for the total pre-equilibrium charm number, 



(he) 2 A 8 / 3 J a 2/ 2 



N M in ~ 4?r ^ 2 ^ y dp ±1 g(p ±1 ) J dp ±2 g(p ±2 ) J dfa j d<p 2 J sa(s) 

! (2tt) 2 [sct(s)] ~ 0.10 (41) 



Therefore NMin/Nsj ~ 1. From the above estimate we can see that although the minimally- 
correlated case allows particles with different rapidities to interact, the dominant contri- 
bution still comes from the region where the two gluons have almost the same rapidity, 
thus there is no sizeable enhancement in the pre-equilibrium charm yield. The Minimally- 
Correlated case is very much like the Bjorken case in that the dominant contribution comes 
from yi ~ y 2 region. 

As a comparison to the above rough estimates in this section, the numerical integration 
gives N Fb = 3.8, N Bj = 0.093, and N Min = 0.078, so N Min /N Bj ~ 80%. 
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IV. DISCUSSION AND SUMMARY 



In this paper, we calculated initial and pre-equilibrium charm production in nuclear 
collisions to test the sensitivity of this probe to the unknown initial conditions in such reac- 
tions. For the initial charm production, the sensitive dependence on the choice of structure 
functions, the Q 2 scale, and the K-factor was noted. The parameters were fixed by fitting 
the limited available experimental data at lower energies. We emphasized the need for new 
measurements of pp and pA charm production to reduce the present large theoretical un- 
certainties. We argued that the copious charm production predicted in ref. || was mainly 
due to the neglect of the coherent suppression of flavor excitation processes. Our calculated 
initial charm yields are close to those computed in ref. || and to the curve 2 in Fig. 2 from 
ref. i. 

For the contribution from pre-equilibrium charm production, we studied the effect of 
correlations between the rapidity y and space-time rapidity rj of mini-jet gluons. For the 
ideal Bjorken-correlated case, where rj = y\ = y2, the pre-equilibrium charm production is 
negligible compared with the yield due to initial gluon fusion. For the opposite extreme 
fireball case, corresponding to uncorrelated y and r], the pre-equilibrium charm production 
is almost a factor of 50 larger than in the Bjorken-correlated case and is comparable with the 
initial charm yield 0. By the estimates of the total pre-equilibrium charm number, we found 
the the difference mainly comes from the r] — y correlation. Therefore, the pre-equilibrium 
charm production is very sensitive to the (rj — y) correlations in the initial state. 

In order to investigate the effect of more realistic correlations that may exist in the initial 
mini-jet plasma, we introduced a minimal correlation model taking into account the uncer- 
tainty principle along the lines of ref. §J. Our main result is that this minimal correlation is 
similar to the ideal Bjorken correlation case and produces negligible pre-equilibrium charm 
compared with the initial charm yield. We also found that the pre-equilibrium charm yield 
is rather insensitive to the formation physics because the early-formed p± > 1 GeV gluons 
dominate. 
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FIGURE CAPTIONS 



Fig.l The cross section for pp — > ccX is plotted as a function of Pi a i>. The solid line is our 
result with M c = 1.3 GeV, K = 3, Q 2 choice-1 and GRV-HO set. The long-dashed 
curve is the result with the same parameters as in Fig. 1 of [23], but using a K-factor 
of 3 instead of doing 0(a 2 s ) calculation. 

Fig. 2 (dN c5 /dY) Y=Q , rapidity density of charm and anticharm pairs for Au — Au collisions 
vs y/S/A. Curves 1-4 from the calculation of Parton Cascade model are compared 
to our calculation of the yield (the solid line) due to initial fusion processes. The top 
curve4 is the total charm production with QGP formation including the incoherent 
flavor excitation processes. Curve3 shows the charm production in the case of QGP 
formation without excitation processes. The bottom curvel is the parton model result 
extrapolated to AA from pp using the A 0J6 scaling measured at much lower energies. 
Curve2 is the parton model result scaled by A 4 / 3 . Our curve uses the asymptotic 
A 4 / 3 scaling. As shown by the two arrows, curve4 becomes curve3 when the coherent 
cancellation of flavor excitation processes is considered, and curvel becomes curve2 
when the high energy scaling is used. So the net dynamical enhancement in the PCM 
(by comparing curve3 to curve2) is comparable to the result of ref. [0]. 

Fig. 3 The production of charmed hadrons as a function of Xf for p—p collisions at Pi a b = 400 



GeV |23]. The solid curve is our result for da/dxj using the first parameterization. The 
dashed curve is our result using the second parameterization. These curves assume a 
delta function charm fragmentation function. 

Fig. 4 The cross section for pp — > bb + X vs y/S/A. The data point at \/S = 630 GeV is 
from ref. | 3"4]| . The dashed cross at =1.8 Tev is obtained indirectly from [23j, and 



the error bar is only illustrative. 
Fig. 5 The mini-jet gluon distribution A~ 4 / 3 (dN/dydp~±) is plotted. The solid curve is 
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taken from the HIJING calculation with radiation effects included, and the circles are 
our result from the initial production. The dashed line is the fit 0.06e~ L25p± . 

Fig. 6 The distribution (Ed 3 N/d 3 p) y=0 of charm quark production using 5(rj — y)-correlation 
is plotted as a function of p±. The solid curve is the initial charm production. The 
curve labelled with filled diamonds is the pre-equilibrium contribution including both 
the soft (p± < 2 GeV) and hard (p± > 2 GeV) components of the mini-jet gluons. 
The curve labelled with unfilled diamonds is the pre-equilibrium contribution including 
only the hard component. 

Fig. 7 The distribution (Ed 3 N/d 3 p) y=0 of charm quark production for the Uncorrelated case 
is plotted as a function of p±. The solid curve is the initial charm production. The 
curve labelled with filled circles is the pre-equilibrium contribution including both the 
soft (p± < 2 GeV) and hard (p± > 2 GeV) components. The curve labelled with 
unfilled circles is the pre-equilibrium contribution including only the hard component. 

Fig. 8 The energy density at z = is plotted as a function of proper time assuming minimal 
correlations and Lorentzian formation probability. The solid curve includes both soft 
and hard components while the dashed curve is calculated using the hard distribution 
and includes only the hard component. 

Fig. 9 The distribution (Ed 3 N/d 3 p) y=0 of charm quark production using minimal r] — y 
correlations is plotted as a function of p±. The curve labelled with filled squares include 
both components while that labelled with unfilled squares include only the fusion of 
hard gluons. The curve labelled with diamonds shows the contribution from fusion of 
soft gluons both with p± < 2GeV. This shows that the pre-equilibrium contribution 
mainly comes from the fusion of soft and hard gluons. 

Fig. 10 The distribution (Ed 3 N/d 3 p) y=0 of charm quark production using different 
formation-time probability distributions. The solid curve is obtained using the 
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Lorentzian form in equation fll8|), and the dashed curve using the theta function form 
in equation (|20|). 
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